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i 
ABSTRACT 


In Calculus we learned that S>y_, k = na). fhat- se ake = nent Gnt)) and 
that.) ak? = (niet) )2, These formulas are useful when solving for the area below a 
quadratic or cubic function over an interval [a,b]. This tedious process, solving for areas 
under a quadratic or a cubic, served as motivation for the introduction of Riemman inte- 
egrals. For the over zealous math student, these steps were replaced by a simpler method 
of evaluating antiderivatives at the endpoints a and b. From my recollection, a former 
instructor informed us do the value of memorizing these formulas. At the very least, 
we needed to know that the sum of the first n integers is a polynomial of degree 2 with 
the leading term sn’, the sum of the squares of the first n integers is a polynomial of 
degree 3 with the leading term gn; lastly, the sum of the cubes of the first n integers is 

1 


a polynomial of degree 4 with the leading term ani. This resulted in the recognition of a 


pattern; hence, the generalization; the sum of the r-th powers of the first n integers is a 


1 


me Thus, we arrive at the theorem 


polynomial of degree r+1 with the leading term 


Theorem 0.0.1. Let r > 0 be an integer. Then there is a unique polynomial, p,(x) of 


degree r +1 with the leading coefficient HI such that for every integer n, 


S> kp): 
k=1 


To some math scholars, these polynomials are called Faulhaber polynomials, 
named after Faulhaber, a German mathematician, who was one of the first mathemati- 
cians to recognize that the sum of r-th powers is indeed a polynomial. His discoveries 
resulted in “simple” forms of formulating these polynomials when r is odd. 

I will prove Theorem 0.0.1, and many other properties of Faulhaber polynomials. 


For example, I will prove that for all r the coefficient of x” is 5. Thus 


_ 1 r+1 i r 
Pr(t) = Hgts ees 


The remaining coefficients of this polynomial, do not have a simple form. This 
is especially the case for the coefficient of x. It turns out that the coefficients of x in 
p(x), are the so called Bernoulli numbers. I will establish this as well as many related 


properties about Bernoulli numbers. 


I will end my thesis with an interesting observation discussed in [Fej05]. The 


formula for the sum of cubes is especially interesting giving us the identity 


13423433 4-0-4 (n—-2)3+(n-18 +n? = (14243454 (n-2) + (n-1) +n)’. 


If in both left hand and right-hand sides of this identity we replace the term (n — 1) by 
2 one would think that the resulting equality 


13423433 +.--+(n—2)? +(2)3 +n? = (14243+4---+(n—2)4+(2) +n) 


is false. It was shown in [Fej05] it is actually true. Moreover the pair ((n — 2), 2) is the 


only such switch that works for all n. 
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Chapter 1 


The Sums 


1.1 Carl Friedrich Gauss and the sum of the first n integers 


Carl Friedrich Gauss, (1777-1855), born in Brunswick Germany, was a child 
prodigy considered to be the greatest German mathematician of the 19*” century [Wel19]. 
He was dubbed the “Prince of Mathematics” and his work influenced a great deal of the 
math areas studied in this day and age which include: number theory, algebra, differential 
geometry among other fields. At age three, he corrected an error in his father’s payroll 
calculations. By age 7, his arithmetic abilities had grown so much so that his instructors 
admitted that there was nothing more that they could teach him [Wel19]. 

In one of his classes, he was given what was meant to be busy work. The prob- 
lem: find the sum of all the numbers from 1 to 100. He solved this in a matter of minutes. 
How? Obviously, individually adding the numbers 1 to 100 is incredibly time consuming. 


If we start by adding the smaller numbers, it’ll look something like this: 


14243 = 6 
14+24+3+44 = 10 
1424+34+4+5 =15 
1424+34+4+4+5+4+6 = 21 
1424+34+4+54+6+7 = 28 
1424+34+4+4+54+64+7+8 = 36 


and so on [Wel19]. 
We don’t know how Gauss calculated the sum of the first 100 numbers; but, he 
may have recognized that when there are even number of numbers to add as in 1 + 2+ 


3+4+54+6+7+8, we could pair 1 with 8, 2 with 7, 3 with 6, and 4 with 5. Each pair 


adds up to 9, and there are four pairs so that the sum is 9-4 = 36. Similarly, if we want 
to add the first 100 numbers, we could pair 1 with 100, 2 with 99, 3 with 97, and so on 
with the last pair being 50 with 51. Each pair sums to 101 and there are 50 pairs. So, 
the sum of the first 100 numbers is 101 - 50 = 5050. 

Another way to add the first 100 integers, with the little help of symbolic algebra, 


is as follows. Let S' denote the sum. Then, 


S= 1 42 43 +---+ 98 +499 +4100 
S= 100 +99 498 +--+ 3 42 41 
2S = 101 Ol, 01 se AO OL ahd 


The first row is just the definition of S. The right hand side of the second row is 
the sum of the first 100 integers written backwards, so it is also S. Finally the third row 
is obtained by adding the corresponding numbers of the two rows above. The left hand 
side of the third row is 2S, while the right hand side is 101-100. Hence 2S = 101 - 100, 
and from here S = tO 100 = 5050. 


This approach can be easily modified to derive the formula for the sum of the 


first n integers. Let S denote the sum. Then, 


S=1 +2 +3 +---+ (n—2) +(n-1) +n 
Si +(n—-1) +(n—-2) +---+ 38 +2 +1 
25=(n+1) 4+(n4+1) 4+(m41) +--+ (n¥#1) +(n+1) +(n+1) 


Similar to that above, the first row is just the definition of S. The right hand 
side of the second row is the sum of the first n integers written backwards, so it is also S. 
Finally the third row is obtained by adding the corresponding numbers of the two rows 


above. The left hand side of the third row is 2S, while the right hand side is (n+ 1)-n. 


Hence 2S = (n+ 1)-n, and from here S = ae 


This proves the following theorem. 


(n+1)-n ; 


Theorem 1.1.1. The sum of the first n integers is —3 


We could generalize this approach to derive the sum of the first n terms of the 


arithmetic progression ao, a9 + d, ag + 2d,...,a9 +(n—1)d. Let S denote this sum. Then 


S =a + (ao +d) (ag + 2d) +++ +(ao + (n—1)d) 
SSG ReaD. ap De) Ape Ga aidi- hoe hae, 
25 =(2a9 + (n—-1)d) +(2a9+(n—1)d) +(2a9+(n—-—1)d) +--+ 4+ (2a9+ (n—-1)d) 


Again, the first row is just the definition of S. The right hand side of the second 
row is the sum of the first n terms written backwards, so it is also S. Finally the third 
row is obtained by adding the corresponding numbers of the two rows above. The left 
hand side of the third row is 2S, while the right hand side is (2a9 + (n — 1)d)-n. Hence 
2S = (2a9 + (n — 1)d)-n, and from here S = Cases, 


Thus we arrive to the following result. 


Theorem 1.1.2. The sum of the first n terms of the arithmetic progression ag,ao + 


.. (2a0-+(n—1)d): 
d,aj + 2d,... is Se 


The special case of Theorem 1.1.2 with ag = 1 and d = 2 is an interesting result 


in its own. 


Corollary 1.1.3. The sum of the first n odd numbers is n?. 


1.2 The sum of squares of the first n integers 


For the sum of the first n squares, the pattern does not hold. Let 


Sp =12 422 4382 $.04 ne? 


=1 +44 49 4---4+ n?. 
In order to obtain the formula for S2, we may be tempted to repeat the same procedure 


that was successful in the aforementioned computation of: 
Sy =14+243+---4+n. 


Thus, write the sum twice, reversing the order the second time: 


So =1 +4 +9 +--+ (n-2)? 4(n-1)?) +n? 


Soar" $(m=1? +(n-2? ++ 9g +4 +1 


Adding both sums results in the sum, 


299 = [1+ n7] + [4+ (n—1)7] + [(m— 2)? +9] +--+ [(n— 1)? +4] + [n? +] 


For this method to work in the right-hand side of this equation all the expressions 
inside the brackets should be identical which obviously is not the case. Thus, we have to 
consider an alternative approach. First we will redo the formula for the sum, 5S}, of the 


first n integers. We begin by writing 


(wel? Sn? + 2n +1 
oe Oe +2(n—1) +1 
(w=1P = (n=27 + 2(n — 2) +1 
5? = 4? +2-4 +1 
A? = 3? +2-3 +1 
3? = 9? ete?) +1 


= 1? Ae A 


After adding these n equations and canceling the numbers 2?,37,...,n? that 


show up to the left and to the right of “=” sign, the resulting equation becomes 
(n +1)? =174+2(1424---+n)+1-n 
or 


(QPS T8-6i sb 


Solving for S; we get 


(n+1)?-—n-1 


Si = 5 


Mimicking this approach for S2, the sum of the squares of the first n integers, 


we use the identity (a + 1)? = a? + 3a? + 3a + 1 to write 


(n+1) =n? + 3n? + 3n +1 
n® = (n— 1)° 3a)? +3(n—1) +1 
(n — 1)? = (n—2)° + 3(n — 2)? + 3(n — 2) +1 
a3" +3-3? +3°3 44 
coe +3-2? +3-2 +1 
2° = 18 A531" +3-1 Al 
After adding these n equations and canceling the numbers 2?,3°,...,n? that show up to 


the left and to the right of “=” sign, the resulting equation becomes 
(n+1PP =P +304 +3 pak +1-n 
or 
(n+1)? =14+35)4 35, +n. 


Using that S; = te, and solving for Sz results in 


_ (nt1)8-1-3 2) nn 
S52 = 3 


which simplifies to the familiar formula 


Sy = n(n + 1)(2n+ 1) 
6 


Thus we have prove the following theorem 


Theorem 1.2.1. The sum, S29, of squares of the first n integers is Sg = mint Y Grr) 


1.3. Sums of cubes and fourth powers of the first n integers 


In this section we will find the formulas for S3, the sum of cubes of the first 
n integers, and for $4, the sum of fourth powers of the first n integers. The method is 
identical with obvious changes to the method used in the previous section to find the 
formula for Sb». 

In order to find the formula for $3 we begin with the identity (a + 1)4 = a4 + 
4da® + 6a? + 4a + 1 to write the following equations 


(n+1)4 =n4 + 4n3 + 6n? +4n +1 
n* =(n—1)* Ain — 1)? + 6(n — 1)? + 4(n — 1) =e | 
(n —1)* = (n—2)* + 4(n — 2)3 + 6(n — 2)? + 4(n — 2) +1 
Cs +4-39 +6-3? +4-3 +1 
ee ides 2F +6-2? +4-2 ae 
geal +4-13 +6-1? +4-1 +1 
After adding these n equations and canceling the numbers 24, 3*,...,n* that show up to 


the left and to the right of ”=” sign, the resulting equation becomes 
(Qt) 1S st ak PO a eee ele 
or 
(n+ 1)* =1+453 +652 +451 +n. 


Using that S; = oy So = oe) and solving for $3 results in 


+1)4 1 gretD@n+)) 4nintd) 


sg = 6 


n 


which can be simplified to 


n?(n +1)? 
a 


Thus we have proved the following theorem 


S3 = 


2 
Theorem 1.3.1. The sum, S3, of cubes of the first n integers is S3 = [2] . 


In order to find the formula for $4 we begin with the identity (a+1)° = a°+5a4+10a? + 
10a? + 5a+ 1 to write the following equations 


(n+1)? =n° + 5n4 + 10n? + 10n? + 5n +1 
Ww =(n—1) 5(n—1)* +10(n— 1)? 10ST 5G = 1D) Ea 
(n — 1)® = (n—-2)° 5(n —2)*  +10(n — 2)3 10(n—2)7 +5(n—2) +41 
4° = 3° aN Se 033° +10-3? +5-3 +1 
3° = 2° eo +1082? Ai" 45:2 +1 
2° =1° eee Te" se Tosie +5-1 za 
After adding these n equations and canceling the numbers 2°,3°,...,n° that 
show up to the left and to the right of ”=” sign, the resulting equation becomes 


(n +1)? =1°+5 7, k*4+ 10 7, P+ lO, P45 Vk +10 
or 


(n+1)° =1+594 +1053 + 10S. + 5S, +n. 


2 
Using that S; = n(n) So = mint Gnt and $3 = ead to solve for $4 results in 


(n+1)® of newt) |” 19 ™tVEn+1) 5 mint) Pee | 


S4= : 


which can be simplified to 


n(n + 1)(2n + 1)(3n? + 3n — 1) 
= ag 


Thus we have proved the following theorem 


Theorem 1.3.2. The sum, S4, of the fourth powers of the first n integers is Sq = 
n(n+1)(2n+1)(3n?+3n—1) 
30 


1.4 Sum of the r-th powers of the first n integers is a poly- 


nomial 


The derivation of formulas for Sg, S3 and S4 followed the same pattern, so it 
is no surprise that these methods can be used to find a formula for S;,, the sum of r-th 


powers of the first n integers. 


Theorem 1.4.1. Let r > 1 be an integer. For1<k<_r, let S, denote the sum of k-th 
powers of the first n integers. Then 
GH ate ye oe 
r+1 ; 
Proof. We begin with the binomial identity (a + 1)"*1 = a™t! + Cea + Cars + 


S, = 


sep (a + 1 to write the following equations 


es OO Oe cas (Ens 
memo +(T man +(e )@-yrt tet (Pen 41 
(m-ytem-yt + (FP )m-ar + (FP )@-ayr (T’)@-2 41 
nO) ae) ce) er 
Se 0) 10) ae oe ae 
ee Ce Ce ae 
After adding these n equations and canceling the numbers 2”+!, 3"+1,...,n"*+1 that show 


up to the left and to the right of “=” sign, the resulting equation becomes 


pel. r+1\c *: r+i\c yo] r+i\c 
(n+1) =1+/ >>: +( cap S ie os (eee Pp Bee 
k=1 k=1 k=1 
or 


1 1 
mit ar4(r4ns+ ("3 Jsatee( )si+n. 


Solving the last equation for S; proves the theorem. 


Corollary 1.4.2. S;,. is a polynomial in n of degree r+1 with the leading coefficient =: 


Proof. The proof is by induction on r. For r = 1, we have seen that S; = sn? + sn 
so the induction case r = 1 is established. Now assume that r > 1 and that for all 
1<k<r-—1, S; is a polynomial of degree k +1. Since by induction hypothesis every 
term in the right-hand side of the formula for S; in Theorem 1.4.1 is a polynomial in n, 
with the term (n+ 1)"*! being of the highest degree, S$, is also a polynomial of the same 

1 


degree as (n + 1)"+! and its leading coefficient is —-. 
r+l 


10 


Chapter 2 


The Numbers 


2.1 Bernoulli numbers 


In this chapter we will take a different approach. We will not assume any of 
the results from Chapter 1. Let r > 0 be an integer and let P(r) denote the set off all 
polynomials p(x) with p(0) = 0 that satisfy the identity 

p(a +1) — p(x) = (a +1)’. (2.1) 


We will show that P(r) has exactly one element which we denote by p,(x). Moreover, we 


will show that p,(n) = S, = \¢_, k” for alln EN. 


Lemma 2.1.1. Let r > 0 be an integer. Suppose that p(x) € P(r). Then there is a 


constant A such that 
q(x) = (r +1) [vat + Ax € P(r+1). 
Proof. Let A be a constant and let 
fey Gad) 7 "BbaE ae, 


Since by assumption p(x) is a polynomial, and since polynomials are closed under inte- 


gration, q(x) is also a polynomial. Note that q(0) = 0. The derivative of q satisfies 


d(x) =(r+1)p(x) + A. (2.2) 


11 


Hence 


(q(@ +1) — g(a)!’ = (+1) -d (2) 


= (r+ 1)(p(x + 1) — p(x) by (2.2) 
=(r4+1)(e4+1)’. by (2.1). 
By taking anti-derivatives we obtain 
q(a@+1)-—q(z)=(@+1)+C (2.3) 


for some constant C. If we take A = 1 — (r+ 1) i p(t)dt then q(1) = 1. Substituting 
x = 0 into (2.3) we obtain 
q(1) — q(0) =1+¢C 
1-0=14+C 


It follows that C = 0. Thus for this choice of A, g(x) € P(r +1). 


Theorem 2.1.2. For every non-negative integer r, P(r) has exactly one element, p,(x). 


Moreover pr(n) = Sop, k” for alin EN. 


Proof. The proof that P(r) 4 @ is by induction on r. For r = 0, it is straight forward to 
check that 2 € P(0). Now assume that for some r > 0, P(r) 40. Let p(x) € P(r). Then 
by lemma 2.1.1, P(r +1) 4 @ and the induction is complete. 

We also use induction this time on n to show that if p(x) € P(r) then p(n) = 
ype k" for all n € N. For n = 1 we have to show that p(1) = 1. Indeed, from (2.1) 
p(0+ 1) —p(0) = (0+1)’, and since p(0) = 0 we get p(1) = 1. Now assume that for some 
n> 1, p(n) = R_, kr. From p(n +1) = p(n) + (n 4:1)" = Ry e+ (n 4.) = OT ke 
completing the induction proof. 


Now if t(a) € P(r) then p(n) — t(n) = 0 for all n € N and since p(x) — t(x) isa 


—t 
polynomial this is possible only if p(a) — t(2) = 0 for all « € R. Hence P(r) has only one 


element. 


From lemma 2.1.1 we have that p,41(#) = (r+1) fj p(t)dt+Az for some constant 
A. From (2.2) we see that A = p/.,,(0). These constants, {A, = pj.(0)}?29 are known 
in literature as Bernoulli numbers. We will show that Bernoulli numbers, {A,}, can be 
generated independent of polynomials, {p,(x)}, and that coefficients of these polynomials 


can be expressed in terms of Bernoulli numbers. 


12 


Lemma 2.1.3. For k > 2, the k-th derivative of pl, (a )=(r+1)pr es D(x). 


Proof. This is immediate from (2.2). 


In the statement of next lemma for k > r+1, the subscript of A;+1_,% is negative, 
and Bernoulli numbers are not defined for negative subscripts. We could restrict the value 
of k tok <r+1 to avoid this, or we could simply set A; = 0 for negative t, and that is 


what we are going to do. 


Lemma 2.1.4. Let k > 1 be an integer. The k-th derivative of p,(x) at 0 satisfies the 
following formula 


kk ort 


p*)(0) 1 /f/r+l1 
k Arti1—k- 


In particular pe Po) =0 fork >r+l. 


Proof. For k = 1 this formula just means that p/.(0) = (TT) Arpt = A, which is 
true by definition, so it needs only to be verified for all & > 2. The proof is by induction 
on r. For r = 0, po(x) = x, so that all its derivatives but the first are zero. We also have 
that for k > 2 Aopii—~ = 0 by definition and the formula is ase for r = 0. 

Now assume that for some r > 0 and for all k > 1, pe Oe a Ape 


r+1\ k 
Let k > 2. By lemma 2.1.3 


po), (0) r+1p*-) 0) 


m = k (ko 0) and by induction hypothesis we get 
r+1 1 r+1 
= es ed oss Apt ied ie) 


1l/fr+l1 1 r+2 
" =; = _ Artank = —s( k Aro 


and the induction is complete. That bE pe *(0) J=Ofork>r+ 1, follows from the formula 


and the definition of A; for t negative. 


In the proof of Theorem 2.1.5 below we will use the fact that po(x) = x, and 
pi(x) = 527 + $2 so that Ap = 1 and Ay = 3 


Theorem 2.1.5. The polynomials p,(x) can be expressed in terms of Bernoulli numbers 
as p,(x) = rar oie (TNA rti-k «*. In addition the two leading terms of p,(x) are 


1 rti af 
ri and 5 gv 
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Proof. Since p,(0) = 0, by Taylor formula 


d {pr (k 
egree of p ph (0) ‘ 


p (2) = > “Fi ae 


k=1 


pi" (0) pet) (0) _ 4 1 
By lemma 2.1.4, +> = 0 for k > r+1, and a Crs ial 7qi4o = 747 80 that degree 
(k) 
of p, isr +1. Also by Lemma 2.1.4, Pr (0) — to (G)) Artie proving the first part of 


Theorem. 


(r) 
For the second part of the Theorem it just remains to show that Pr (0) = 5: To 
that end by Lemma 2.1.4 with k = r we get 


(r) 
pr (0) _ 1 r+l = = 
rl al r ) Artic =A 


1 
a: 


Theorem 2.1.5 and the fact that p,(1) = 1 for all r, allows us to express Bernoulli 


numbers recursively independent of p,(x). 


Corollary 2.1.6. Bernoulli numbers can be defined recursively as Ag = 1 and A, = 


1 r—-1 (r+l 
ea aa : )An forr>1. 


Proof. From p,(1) = 1 and Theorem 2.1.5 we have 


r+l 
1 r+1 
1 =p,(1) = pe 
m= ("gare 


using substitution n = r+ 1—k, and that (2 = (aed) = ce), we can rewrite this 


as 


Solving for A, finishes the proof. 
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Here is the list of the first twenty one Bernoulli numbers obtained using Math- 
ematica code: 
Input: 
A{0]=1 
Table[A[r] = 1 - 1/(r + 1)*Sum|[Binomial[r + 1, n]*A[n], {n, 0, r - 1}], {r, 0, 20}] 
Output: 
1, 1/2, 1/6, 0, -(1/30), 0, 1/42, 0, -(1/30), 0, 5/66, 0, -(691/ 2730), 0, 7/6, 0, -(3617/510), 
0, 43867/798, 0, -(174611/330) 
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Here is the list of p,(x) for r = 0,1,2,...,20. 


Po xr)=—=z 
i 1 
v)=-2# 4+ =a 
Pl 3 2 
1. 1 1 
228 ee) Line a 
p2(x2) = au 3% rd 
_l4 1 3 1 5 
p3(x) = a rl a 
ee t.§ iv 1 
= + =a 4+ =x -—2£ 
UE ie 2 3 30 
1 1 5 1 
= OS 6 a) Mf eo pe gd To ae 
BE ge 2” 12” 12 
1s Lag Lo 1% 1 
— aes + —-— + —-— ae + —2x 
cis Sa 2” 2” e 22 
= 1 8 1 e 7 6 ig 4 1 2 
BEE ge 2” 12” 24” + 79* 
it ag ie ahs ate 3 1 
=- + 22 ee = + =x Sa 
BEM og 2 3 15 9 30 
1 10 1 9 308 7 6 lia 3 5 
= — aie =, ab. “a Be: —-—2 
BN ATS 2” 4” 10° +5 20 
lou 1 io 5 9 7 1 3 5 
= — + — + — - +2° —=2£ + — 
SO 2” as a 2 66 
vies Pans hs 6 Li wee ily ou 
Pa - 12 8. hae 8 12 
1.45 ae 29 44 1% 22. B35 5 3 
= + == 2 Se 
eae eae 13° 738° 6 7 10 3 
1 1 1 14. 14 65 691 
pis(w) = a4 4 al ee ee eae era +—a* — 2? 
14 84 30 28 20 12 420 
er 1 ral 14. 14 91 - 691 
Pa eee) 74 sa 39 Pe, snes eels 
15 30 150 18 10 6 90 
1 277 429 455 691 35 
re ee 10 28 78 _ SPL a Pee 
16 672 16 12 24 4 
1 1 1 2 1382 - 140 3617 
Pi6(x@) = gl? att . 2 3 ae — —2 f. gs - ——2 
17 357 255 3 15 3 510 
Hae Lee Lary to 17 10 105 11747 6 | 595 4 
18 12 252 270 6 45 3 
1 23494 z 3617 
pis(x) = x9 f git g gil _ ue _ = —_" + 714x° 2 
19 114 95 627 35 10 
1 1 1 13437877 68723 
Do a PO) ag 89 ci0 = TOARTONT (8° a opi 98 4 
20 40 120 88 54600 40 
1 i 2 68723 - 219335 174611 
p20(x@) = —72t 4 Lg l9 | 2 gi! oueee x + 646027 — ——~2> | x x 


21 9 154 28665 10 63 330 
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Notice that for r > 1 and r odd, all A,’s in the list are zero. That this is always 


true is a part of Theorem 2.1.7 below. 


Theorem 2.1.7. For r odd, we have p,(x) — p-(—x) = x". It follows that A, = 0 for 
r > 1 and odd. 


Proof. First we will show that for r odd, p,(x) — p,-(—x) = x” for all x € R. Since 
p,(x) — pr(—2) is a polynomial it is enough to verify this for « = n for alln € N. We 
begin by substituting « = n—1, x =n—2,--- ,x = —n into (2.1) to obtain the following 


identities: 


Pr(—(n — 1)) = pr(—n) — (n — 1)’. 


In the bottom half of these identities we used that r is odd so that (—k)” = —k". By 
adding these identities, and taking into account that all the terms except p,(n), n” and 
pr(—n) cancel, we obtain p,(n) = p-(—n) +n” or p(n) — pr(—n) =n". 

Next we will show that A, = 0 for r > 1 and odd. ‘To that end we first 
differentiate p,(x) — p,(—x) = x”, to get p).(x)+p).(—x) = ra’. Substitution z = 0 into 


the last equation yield A, = p/.(0) = 0 for r > 1 and the proof is complete. 


Corollary 2.1.8. [fr > 1 is odd, then p,—1(x) = Pr(a) | 


ie 
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Proof. Let r > 1 and r odd. Let p(x) = Pr(@) In order to show that p,—1(x) = p(x) 


we need to show that p(x) € P(r —1). That p(x) is a polynomial is obvious. By 
Theorem 2.1.7, A, = 0 so that p(0) = 0. It remains to show that p(x + 1) — p(x) = 
(x +1)"—!, but this follows immediately by differentiating p,(a2 +1) —p,(#) = (x +1)" 


and dividing the resulting equation through by r. 


Corollary 2.1.9. [fr #0 even, then p,(x) + p,(—x) = 2". 


Proof. By Theorem 2.1.7 pp4i(x) — ppyi(—xv) = 2"t!. Now differentiating the last 


equation we get pi.4(x) + pi4,(—2) = (r+ 1)a". Dividing through by r + 1 we ob- 


Prti(t) | Prpi(—®) ar P41 (2) 


ee ea x". But by corollary 2.1.7 p-(7) = —54- proving that 


tain 


Pr(x) + p,(—2) =x". 


We already know that p,(x) is a polynomial of degree r + 1 and that the coef- 


ficients of "+! and 2” are HT and 5 respectively. The next corollary shows that “every 


other” coefficient is zero. 


Corollary 2.1.10. For r odd, p,(x) — aq — $2" is of the form eer aja" 1-2, 


, 1 1 j 2-94 
while for r even, p,(x) — are — 52x" is of the form j=0 ang A= A8 


Proof. For r odd we have to prove that p*)(0) = 0 for k # r and k odd. This follows 
by differentiating k-times the identity p,(x) — p,(—x) = x” from Theorem 2.1.7 to obtain 
pr”’(0) + pr(k)(0) = 0. 

For r even we have to prove that p*)(0) = 0 fork # r and k even. This follows 
by differentiating k-times the identity p,(x) + p-(—x) = x” from corollary 2.1.7 to obtain 
pS” (0) + p (0) = 0. 
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2.2 Faulhaber 


Bernoulli’s formula is sometimes called Faulhaber’s formula. Johann Faulhaber 
of Ulm was a notable German mathematician whose research was particularly significant 
in providing alternative methods of computing the sums of powers. In his findings he 
recognized the patterns in the sums which led to the formulation of polynomial identities 
that express the sum powers of the first n integers. The resulting polynomials have 
been extensively studied by a multitude of math scholars. However, Faulhaber’s main 
contribution to the study of these polynomials was that he was the first one to realize 
that for r odd, the coefficients of p;(x) can be expressed in terms of x(x +1). We will 
end this chapter by proving this. In what follows we are not assuming any results that 
we proved so far. 

First we need the following lemma regarding odd polynomials. A polynomial, 
p(x) is odd if p(—x) = p(x). Since odd polynomials are closed under addition and scalar 
multiplication, the set of all odd polynomials, O, is a subspace of the set of all polynomials. 
The natural basis for the vector space of odd polynomials is {a?/ ocd ee It is also true 
that the set of all odd polynomials of degree less than or equal to 2k — 1 is a subspace of 
O with the dimension equal to k. 


Lemma 2.2.1. Let r be an odd integer. The sequence {x/ ((a + 1)! — (a — 1)’) jai Of 


odd polynomials forms the basis for all odd polynomials of degree less than or equal to r. 
Proof. Let qj(x) = x) ((a + 1)/ — (z — 1)). Then 
qj(—a) = (—a)’ ((-a + 1) — (-a — 1) 
Cin ea T= sary) 
= 2) ((@—1) — («+ 1)’) = -a(2) 


proving that q;(x) are odd polynomials. The degree of q;(x) is 27 — 1, so that the set 
41 


W= {aj}j24 is linearly independent, and since W has rt elements it is the basis for 


r 


the space of all odd polynomials of degree less than or equal to r. 


If r is an odd integer then x” is an odd polynomial so by Lemma 2.2.1, 2” is a 
r+ 
2 


linear combination of {x ((a + 1)? — (x — 1)’) jel: 


Theorem 2.2.2. Let r = 2m—1. If x” = ofl, aga® ((a +1) — (a — 1)*), then p,(x) = 
TM anh (w+ 1). 


Proof. Let p(x) = S77, agt®(a + 1)*. Then p(—x) = S77, ax(—1)¥a*(—1)* (2 — 1)* so 
that p(x) — p(—«) = 
Dore ant* ((@ + 1)* — (a — 1)*) =". In particular 


p(a +1) — p(-(a + 1)) = (@ +1)’. 


But p(—(# + 1)) = Sophy ae(—(@ +: 1))*((-(@ +1) +18 
= 1 an(—1)*(@ + 1)*(-1)¥(x)* = p(x). Hence p(x + 1) — p(x) = (2 +1)". Since 
p(0) = 0, p(x) € P(r), completing the proof. 


Thus for r odd, Faulhaber polynomial can indeed be expressed as a polynomial 


in z(# +1). For r even by Corollary 2.1.8 we have the following result. 


r+2 r 
Corollary 2.2.3. Letr be even. If pp4i (x) = 0,2, anv’ (x+1)*, then p,(x) = aut 70 
1)a} (a +1). 


1 r+2 
Proof. By corollary 2.1.8, p,(a) = ae Sy But piii(z) = £30,2, an(2? + 2)*, so by 


r+2 
Chain rule, p},j(x) = (2a +1)>0,2, axk(z? +. 2)*-1. Let j = k — 1, then pi, ,(z) = 


(2x + 1) ee, aj4i(j +1)(a? + x). Dividing through by r + 1 completes the proof. 
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Here is the list of p,(a) expressed in the terms of x(a +1) for the first 10 odd 
values of r. The list is generated using C++ codes illustrated in Appendix A. 


pila) = 50(0 +1) 

pala) = 4(2(@ +1)? 

pola) = sleet 0))® = 5 (ele +1)? 

m(a)=s(e(e+ 0) Flee +) Flee +)? 

roa) = alae +> Feet Hy tele t HS @e +1)? 

pula) = sleet 1))® — —Z(el@+1)® + Greet ))t — — ale + 1)) 

f 2 o(a +1))? 

pis() = (eet 1) eet 1S +2 ete  — - P@@ +1) 
sept) — + Foa(o(e+ 1)? 

pis() = ze(o(et+1)® — —F(ale +1)? — + Ewe +)® — — See +19) 
Set) -Pe@ety + 2e@e@+1)? 

pile) =e +9) - Feet +2e@+y? = -Be@w@+1)° 
SS ee+y -See@+p +S %@@+ns  -Sle@+y) 

pio(e) = sa(ale +i) Set) + ele tS - ee +1)" 
Sale +i)®  - e+ +2 eo ery - SH ew@+1) 
43867 
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Chapter 3 


The Discovery 


3.1 On sum of cubes 


We have found formulas, S,, for the sums of r-th powers of the first n integers 
for the first twenty values of r. An interesting and well-known observation is that $3 = S$? 
that is 

Pf Phot gS 1 Sob bay. (3.1) 


What may come as a surprise is that if for n > 2, in (3.1), we replace (n — 1) by 2 on the 


both sides of the equations the resulting equation 


POF es (ye Oy 08 hay = (ee Os ke Gy 2) en)? (3.2) 


is still true. This remarkable identity has been studied in [Fej05]. In this paper the 
authors consider a more general question of finding all triples (k,z,n) with 1 <k,a<n 


for which if & in (3.1) is replaced by « # k the resulting identity 


k-1 n k-1 n 
SOP+e+ SO KP =Ols+e4+ SO GP? (3.3) 
j=l j=k+1 j=l j=k+l 


still holds. The main result in [Fej05] is the following theorem. 


Theorem 3.1.1. Let k, x, and n be positive integers such that] <k,a<n,n4#k. Then 
the triple (k,x,n) satisfies (3.3) if and only if 

(2= Fork) = (ni — 1); 2eny or 

(b) (k,x,n) = (2(p—1) + 3p 2p+s,x+k-—p) where integers p and s satisfy p> 2 
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and s is a positive divisor of 3p(p—1). Moreover, different values of the pair (p, s) yield 


different solutions. 


We can actually use Theorem 3.1.1 to not only find a triple (k,z,n) satisfying 
equation (3.3) but also to determine if a given triple is in fact, a solution. 
For example: if p = 20, then we calculate 3(p — 1)p = 1140, and for s = 10 we obtain 
the triple (152,50, 182) while if we pick another positive divisor of 1140, say s = 19, we 
obtain the triple (98, 59, 137). 


Hence, 
151 182 151 182 
LP +503 + S > 73 = 273,935, 601 = ( PO a Seay 
jJ=153 j=153 
and 
137 137 
a +593 + S~ 7? = 88, 623, 396 = ( Ly +e04 Soa) 
j=99 j=99 


On the other i (97, 37, 130) is not a solution because this triple would yield p = 
r+k—n= 98+ 37 — 130 = 5 and thus, s = x — 2p = 37 — 10 = 27 would have to bea 
divisor of 3(p — 1)p = 60. 

When we check the validity of these values we get, 


597, 59 +373 + 3189, 73 = 71,614, 686 
A T1,470, 116 = (2, § +37 + TS, 9). 


Theorem 3.1.1 part b) naturally raises questions about what else can be said 
about the “non-trivial” triples (k,xz,n). Obviously from Theorem 3.1.1 part b) the for- 
mulas for k and x give us the lower bounds for k and z. Namely k > 2(p—1)+1> 
2(2-1)+1=3 and 2=2p+s>2p4+1>5. 

What if k > 3 is given, can we find x and n so that (k,z,n) is a “non-trivial” 
solution to (3.3). Or if > 5 is given, can we find k and n so that (k,x,n) is a “non-trivial” 
solution to (3.3). The answer to both questions is yes. 

First if k > 3 is given, and odd we can take p = it and s = 3(p—1)p so that 
k =2(p—1)+ 3p Dp If k > 3 is given and even we can take p = f and s = 3 ep sO 
that k = 2(p—1)+ alee 

In the case x > 5 is given and odd we can take p = ot and s = 1 so that 


x =2p+s. If x > 5 is given and even we can take p= 5 —1 and s = 2 so that x = 2p+s. 
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An obvious question is what if n is given. Can we find k and x so that (k, x,n) 
is a “non-trivial” solution to (3.3). In this case the answer is yes if and only if n?++n+1 
is not a prime number or 3 times a prime number. |fej05| 

We end this chapter by listing all “non-trivial” solution to (3.3) for all 2< p< 
153. The following table was generated with the help of C++ code detailed in Appendix 
B: 


Table 3.1: Nontrivial Solution Chart 


k x n k x n k x n 
p=2]| 4 4 9 p=8 | 26] 30 | 48 | p=13 |] 33 | 78 | 98 
3) 6 9 22) 37 | 51 28 | 143 | 158 
3 | 10) 11 21) 40 | 53 26 | 260 | 273 
8 5 11 20 | 44 | 56 25 | 494 | 506 
p=3 | 6 | 15 | 18 18 | 58 | 68 | p=14 | 47 | 54 | 87 
13 | 8 18 17} 72 | 81 40 | 67 | 98 
7 | 12 | 16 16 | 100 | 108 39 | 70 | 95 
10; 9 16 15 | 184 | 191 33 | 106 | 125 
6 | 15 | 18 | p=9 | 24} 45 | 60 32 | 119 | 137 
5 | 24 | 26 20) 72 | 838 29 | 210 | 225 
p=4 | 12) 14 9 18 | 126 | 135 28 | 301 | 315 
6 | 15 | 18 17 | 234 | 242 27 | 574 | 587 
10 | 17 | 23 | p=10 | 28 | 47 | 65 | p=15 | 46] 65 | 96 
8 | 26 | 30 23 | 74 | 87 42 | 75 | 102 
7 | 44 | 47 20 | 155 | 165 38 | 93 | 116 
p=5 | 12) 20 | 29 19 | 290 | 299 37 | 100 | 122 
14) 20 | 29 | p=11 | 35 | 44 | 68 35 | 120 | 140 
18} 16 | 29 31 | 52 | 72 33 | 156 | 174 
13 | 22 | 30 30] 55 | 74 30 | 345 | 360 
12 | 25 | 32 26 | 77 | 92 29 | 660 | 674 
11 | 30 | 36 25 | 88 | 102 | p=16 | 50] 68 | 102 
10 | 40 | 45 23 | 132 | 144 48 | 72 | 104 
9 | 70 | 74 22 | 187 | 198 46 | 77 | 107 
p=6 | 10) 22 | 35 21 | 352 | 362 40 | 104 | 128 
20} 21 | 35 | p=12 | 40] 46 | 74 39 | 112 | 135 
15 | 30 | 39 33 | 60 | 81 38 | 122 | 144 
12) 57 | 68 31 | 68 | 87 35 | 176 | 195 
11 | 102 | 107 26 | 123 | 137 34 | 212 | 230 
p=7 | 21 |) 28 | 42 24 | 222 | 234 32 | 392 | 408 
19 | 32 | 44 23 | 420 | 431 31 | 752 | 767 
14) 77 | 84 | p=13 | 42 | 52 | 81 | p=17 | 56] 68 | 107 
13 | 140 | 146 37 | 62 | 86 49 | 82 | 114 
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Table 3.2: Nontrivial Solution Chart 


k x n k x n k x n 
p=17 | 48) 85 | 116 | p=22 | 64 | 107 | 149 | p=26 | 60 | 247 | 281 
44 | 102 | 129 60 | 121 | 159 56 | 377 | 407 
40 | 136 | 159 56 | 143 | 177 55 | 442 | 471 
38 | 170 | 191 53 | 170 | 201 53 | 702 | 729 
36 | 238 | 257 51 | 198 | 227 | p=27 | 78 | 135 | 186 
35 | 306 | 324 49 | 242 | 269 65 | 216 | 254 
34 | 442 | 459 44 | 737 | 759 | p=28 | 82 | 137 | 191 
33 | 850 | 866 | p=23 | 77 | 92 | 146 68 | 218 | 258 
p=18 | 61 | 70 | 113 67 | 112 | 156 61 | 380 | 413 
51 | 90 | 123 66 | 115 | 158 58 | 623 | 653 
36 | 495 | 513 55 | 184 | 216 | p=29 | 98 | 116 | 185 
35 | 954 | 971 50 | 299 | 326 85 | 142 | 198 
p=19 | 63 | 76 | 120 AT | 552 | 576 84 | 145 | 200 
55 | 92 | 128 46 | 805 | 828 77 | 174 | 222 
38 | 551 | 570 | p=24 | 82 | 94 | 152 70 | 232 | 273 
p=20 | 68 | 78 | 126 69 | 120 | 165 68 | 261 | 300 
58 | 97 | 135 64 | 140 | 180 63 | 406 | 440 
57 | 100 | 137 55 | 232 | 263 62 | 464 | 497 
53 | 116 | 149 54 | 255 | 285 60 | 667 | 698 
50 | 135 | 165 50 | 462 | 488 59 | 870 | 900 
48 | 154 | 182 48 | 876 | 900 | p=30 | 103 | 118 | 191 
44 | 230 | 254 | p=25 | 88 | 95 | 158 87 | 150 | 207 
43 | 268 | 291 84 | 100 | 159 76 | 205 | 251 
42 | 325 | 347 73 | 122 | 170 68 | 321 | 359 
41 | 420 | 441 68 | 140 | 183 67 | 350 | 387 
40 | 610 | 630 66 | 150 | 191 63 | 582 | 615 
p=21 | 76 | 77 | 132 58 | 230 | 263 | p=31 | 105 | 124 | 198 
75 | 78 | 132 57 | 250 | 282 91 | 152 | 212 
68 | 87 | 134 56 | 275 | 306 78 | 217 | 264 
60 | 105 | 144 53 | 410 | 438 70 | 341 | 380 
58 | 112 | 149 52 | 500 | 527 69 | 372 | 410 
54 | 132 | 165 50 | 950 | 975 65 | 620 | 654 
50 | 168 | 197 | p=26 | 89 | 102 | 165 | p=32 | 110 | 126 | 204 
49 | 182 | 210 80 | 117 | 171 94 | 157 | 219 
AT | 222 | 248 76 | 127 | 177 93 | 160 | 221 
45 | 294 | 318 75 | 130 | 179 86 | 188 | 242 
44 | 357 | 380 65 | 182 | 221 78 | 250 | 296 
42 | 672 | 693 63 | 202 | 239 74 | 312 | 354 
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Table 3.3: Nontrivial Solution Chart 


k x n k x n k x n 
p=32 | 70 | 436 | 474 | p=37 | 126 | 148 | 237 | p=42 | 91 | 658 | 707 
68 | 560 | 596 109 | 182 | 254 89 | 822 | 869 
66 | 808 | 842 99 | 222 | 284 | p=43 | 147 | 172 | 276 
p=33 | 108 | 138 | 213 | p=38 | 131 | 150 | 243 127 | 212 | 296 
100 | 154 | 221 112 | 187 | 261 102 | 387 | 446 
96 | 165 | 228 111 | 190 | 263 98 | 473 | 528 
86 | 210 | 263 93 | 298 | 353 93 | 688 | 738 
82 | 242 | 291 80 | 779 | 821 91 | 860 | 908 
80 | 264 | 311 | p=39 | 114 | 195 | 270 | p=44 | 152 | 174 | 282 
75 | 354 | 396 102 | 249 | 312 130 | 217 | 303 
73 | 418 | 458 95 | 312 | 368 129 | 220 | 305 
72 | 462 | 501 94 | 325 | 380 119 | 260 | 335 
68 | 858 | 893 89 | 420 | 470 108 | 346 | 410 
p=34 | 100 | 167 | 233 85 | 572 | 618 98 | 561 | 615 
88 | 221 | 275 | p=40 | 143 | 152 | 255 97 | 604 | 657 
84 | 255 | 305 130 | 170 | 260 | p=45 | 143 | 198 | 296 
83 | 266 | 315 123 | 184 | 267 142 | 200 | 297 
77 | 374 | 417 118 | 197 | 275 132 | 225 | 312 
75 | 442 | 483 114 | 210 | 284 115 | 310 | 380 
p=35 | 119 | 140 | 224 104 | 260 | 324 110 | 360 | 425 
110 | 155 | 230 98 | 314 | 372 108 | 387 | 450 
103 | 172 | 240 96 | 340 | 396 99 | 630 | 684 
102 | 175 | 242 91 | 440 | 491 98 | 684 | 737 
98 | 189 | 252 88 | 548 | 596 | p=46 | 144 | 207 | 305 
89 | 240 | 294 87 | 600 | 647 136 | 227 | 317 
85 | 280 | 330 86 | 665 | 711 117 | 322 | 393 
83 | 308 | 356 | p=41 | 140 | 164 | 263 113 | 362 | 429 
82 | 325 | 372 121 | 202 | 282 100 | 713 | 767 
78 | 427 | 470 120 | 205 | 284 | p=47 | 161 | 188 | 302 
75 | 580 | 620 110 | 246 | 315 139 | 232 | 324 
74 | 665 | 704 104 | 287 | 350 138 | 235 | 326 
73 | 828 | 867 100 | 328 | 387 115 | 376 | 444 
p=36 | 124 | 142 | 230 95 | 410 | 464 | p=48 | 166 | 190 | 308 
105 | 180 | 249 92 | 492 | 543 141 | 240 | 333 
98 | 207 | 269 90 | 574 | 623 130 | 284 | 366 
97 | 212 | 273 88 | 697 | 744 112 | 472 | 536 
90 | 261 | 315 86 | 902 | 947 110 | 519 | 581 
84 | 342 | 390 | p=42 | 145 | 166 | 269 103 | 848 | 903 
80 | 450 | 494 123 | 210 | 291 102 | 942 | 996 
77 | 612 | 653 100 | 371 | 429 | p=49 | 168 | 196 | 315 
75 | 828 | 867 96 | 453 | 507 159 | 210 | 320 


Table 3.4: Nontrivial Solution Chart 


k x n k x n k x n 
p=49 | 152 | 224 | 327 | p=54 | 187 | 214 | 347 | p=59 | 203 | 236 | 380 
145 | 242 | 338 159 | 270 | 375 175 | 292 | 408 
132 | 294 | 377 | p=55 | 189 | 220 | 354 174 | 295 | 410 
124 | 350 | 425 163 | 272 | 380 145 | 472 | 558 
114 | 490 | 555 130 | 515 | 590 | p=60 | 208 | 238 | 386 
112 | 539 | 602 119 | 920 | 984 177 | 300 | 417 
110 | 602 | 663 | p=56 | 198 | 217 | 359 163 | 356 | 459 
105 | 882 | 938 194 | 222 | 360 154 | 415 | 509 
p=50 | 173 | 198 | 321 187 | 232 | 363 138 | 651 | 729 
168 | 205 | 323 180 | 244 | 368 136 | 710 | 786 
148 | 247 | 345 176 | 252 | 372 | p=61 | 210 | 244 | 393 
147 | 250 | 347 170 | 266 | 380 181 | 302 | 422 
140 | 275 | 365 166 | 277 | 387 165 | 366 | 470 
133 | 310 | 393 165 | 280 | 389 156 | 427 | 522 
128 | 345 | 423 154 | 322 | 420 140 | 671 | 750 
123 | 394 | 467 152 | 332 | 428 138 | 732 | 809 
119 | 450 | 519 150 | 343 | 437 | p=62 | 215 | 246 | 399 
113 | 590 | 653 145 | 376 | 465 184 | 307 | 429 
112 | 625 | 687 143 | 392 | 479 183 | 310 | 431 
108 | 835 | 893 140 | 420 | 504 153 | 490 | 581 
p=51 | 185 | 192 | 326 138 | 442 | 524 | p=63 | 186 | 315 | 438 
150 | 255 | 354 134 | 497 | 575 178 | 343 | 458 
145 | 272 | 366 132 | 532 | 608 155 | 504 | 596 
134 | 327 | 410 131 | 552 | 627 151 | 560 | 648 
125 | 408 | 482 130 | 574 | 648 | p=64 | 234 | 240 | 410 
118 | 527 | 594 125 | 728 | 797 190 | 317 | 443 
117 | 552 | 618 124 | 772 | 840 182 | 344 | 462 
110 | 867 | 926 122 | 882 | 948 180 | 352 | 468 
p=52 | 170 | 221 | 339 | p=57 | 188 | 240 | 371 158 | 506 | 600 
154 | 257 | 359 184 | 247 | 374 154 | 560 | 650 
138 | 325 | 411 175 | 266 | 384 153 | 576 | 665 
136 | 338 | 422 168 | 285 | 396 142 | 884 | 962 
128 | 410 | 486 150 | 366 | 459 | p=65 | 232 | 250 | 417 
120 | 546 | 614 148 | 380 | 471 224 | 260 | 419 
119 | 572 | 639 140 | 456 | 539 208 | 286 | 429 
115 | 716 | 779 131 | 618 | 692 206 | 290 | 431 
p=53 | 182 | 212 | 341 130 | 646 | 719 193 | 322 | 450 
157 | 262 | 366 126 | 798 | 867 192 | 325 | 452 
156 | 265 | 368 | p=58 | 172 | 287 | 401 188 | 338 | 461 
143 | 318 | 408 152 | 377 | 471 180 | 370 | 485 
130 | 424 | 501 143 | 458 | 543 176 | 390 | 501 


Table 3.5: Nontrivial Solution Chart 


k x n k x n k x n 
p=65 | 160 | 520 | 615 | p=71 | 245 | 284 | 458 | p=77 | 229 | 382 | 534 
158 | 546 | 639 211 | 352 | 492 228 | 385 | 536 
154 | 610 | 699 210 | 355 | 494 218 | 420 | 561 
152 | 650 | 737 182 | 497 | 608 209 | 462 | 594 
148 | 754 | 837 175 | 568 | 672 196 | 553 | 672 
144 | 910 | 989 170 | 639 | 738 194 | 532 | 689 
p=66 | 240 | 249 | 423 161 | 852 | 942 190 | 616 | 729 
229 | 262 | 425 | p=72 | 250 | 286 | 464 185 | 686 | 794 
220 | 275 | 429 213 | 360 | 501 180 | 781 | 884 
195 | 330 | 459 196 | 428 | 552 | p=78 | 280 | 299 | 501 
185 | 366 | 485 169 | 712 | 809 271 | 310 | 503 
175 | 418 | 527 | p=73 | 252 | 292 | 471 253 | 338 | 513 
156 | 627 | 717 217 | 362 | 506 245 | 354 | 521 
152 | 717 | 803 198 | 438 | 563 231 | 390 | 543 
148 | 847 | 929 171 | 730 | 828 217 | 442 | 581 
p=67 | 231 | 268 | 432 | p=74 | 257 | 294 | 477 180 | 849 | 951 
199 | 332 | 464 220 | 367 | 513 | p=79 | 273 | 316 | 510 
154 | 737 | 824 219 | 370 | 515 235 | 392 | 548 
150 | 871 | 954 183 | 586 | 695 182 | 869 | 972 
p=68 | 236 | 270 | 438 | p=75 | 238 | 335 | 498 | p=80 | 278 | 318 | 516 
202 | 337 | 471 222 | 375 | 522 238 | 397 | 555 
201 | 340 | 473 198 | 483 | 606 237 | 400 | 557 
185 | 404 | 521 193 | 520 | 638 218 | 476 | 614 
168 | 538 | 638 185 | 600 | 710 206 | 555 | 681 
p=69 | 228 | 291 | 450 173 | 816 | 914 198 | 634 | 752 
204 | 345 | 480 | p=76 | 250 | 323 | 497 188 | 792 | 900 
182 | 444 | 557 245 | 332 | 501 | p=81 | 240 | 405 | 564 
172 | 529 | 632 240 | 342 | 506 200 | 648 | 767 
170 | 552 | 653 226 | 377 | 527 | p=82 | 244 | 407 | 569 
159 | 750 | 840 200 | 494 | 618 203 | 650 | 771 
p=70 | 253 | 266 | 449 195 | 532 | 651 | p=83 | 287 | 332 | 536 
228 | 301 | 459 188 | 602 | 714 247 | 412 | 576 
208 | 347 | 485 186 | 627 | 737 246 | 415 | 578 
201 | 370 | 501 175 | 836 | 935 205 | 664 | 786 
184 | 455 | 569 | p=77 | 285 | 286 | 494 | p=84 | 292 | 334 | 542 
183 | 462 | 575 284 | 287 | 494 249 | 420 | 585 
173 | 554 | 657 266 | 308 | 497 229 | 500 | 645 
161 | 770 | 861 236 | 363 | 522 202 | 749 | 867 


Table 3.6: Nontrivial Solution Chart 


k x n k x n k x n 
p=85 | 308 | 323 | 546 | p=91 | 315 | 364 | 588 | p=96 | 262 | 572 | 738 
294 | 340 | 549 310 | 371 | 590 235 | 800 | 939 
287 | 350 | 552 271 | 452 | 632 | p=97 | 336 | 388 | 627 
258 | 408 | 581 250 | 533 | 692 289 | 482 | 674 
253 | 422 | 590 245 | 560 | 714 264 | 582 | 749 
238 | 476 | 629 234 | 637 | 780 | p=98 | 341 | 390 | 633 
236 | 485 | 636 | p=92 | 338 | 345 | 591 292 | 487 | 681 
231 | 510 | 656 320 | 366 | 594 291 | 490 | 683 
213 | 646 | 774 274 | 457 | 639 243 | 778 | 923 
204 | 765 | 884 273 | 460 | 641 | p=99 | 350 | 387 | 638 
203 | 782 | 900 266 | 483 | 657 294 | 495 | 690 
202 | 800 | 917 260 | 506 | 674 273 | 576 | 750 
p=86 | 299 | 342 | 555 251 | 548 | 707 250 | 737 | 888 
272 | 387 | 573 234 | 667 | 809 245 | 792 | 938 
256 | 427 | 597 228 | 730 | 866 | p=100 | 333 | 420 | 653 
255 | 430 | 599 224 | 782 | 914 308 | 470 | 678 
221 | 602 | 737 221 | 828 | 957 306 | 475 | 681 
213 | 682 | 809 | p=93 | 308 | 393 | 608 298 | 497 | 695 
204 | 817 | 935 276 | 465 | 648 253 | 740 | 893 
p=87 | 258 | 435 | 606 246 | 600 | 753 252 | 750 | 902 
230 | 561 | 704 230 | 744 | 881 248 | 794 | 942 
215 | 696 | 824 | p=94 | 280 | 467 | 653 | p=101 | 350 | 404 | 653 
p=88 | 290 | 374 | 576 248 | 611 | 765 301 | 502 | 702 
273 | 408 | 593 233 | 746 | 885 300 | 505 | 704 
262 | 437 | 611 | p=95 | 329 | 380 | 614 275 | 606 | 780 
246 | 495 | 653 302 | 425 | 632 260 | 707 | 866 
232 | 572 | 716 283 | 472 | 660 250 | 808 | 957 
218 | 698 | 828 282 | 475 | 662 | p=102 | 355 | 406 | 659 
210 | 814 | 936 245 | 660 | 810 303 | 510 | 711 
p=89 | 308 | 356 | 575 235 | 760 | 900 | p=103 | 357 | 412 | 666 
265 | 442 | 618 | p=96 | 350 | 363 | 617 307 | 512 | 716 
264 | 445 | 620 342 | 372 | 618 | p=104 | 362 | 414 | 672 
242 | 534 | 687 334 | 382 | 620 310 | 517 | 723 
220 | 712 | 843 285 | 480 | 669 309 | 520 | 725 
p=90 | 313 | 358 | 581 280 | 496 | 680 284 | 620 | 800 
267 | 450 | 627 270 | 534 | 708 258 | 826 | 980 
232 | 625 | 767 266 | 552 | 722 | p=105 | 388 | 392 | 675 


Table 3.7: 


Nontrivial Solution Chart 
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k x n k x n k x n 
p=105 | 348 | 444 | 687 | p=113 | 392 | 452 | 731 | p=122 | 363 | 610 | 851 
338 | 462 | 695 337 | 562 | 786 | p=123 | 366 | 615 | 858 
334 | 470 | 699 336 | 565 | 788 326 | 795 | 998 
325 | 490 | 710 308 | 678 | 873 | p=124 | 410 | 527 | 813 
312 | 525 | 732 | p=114 | 397 | 454 | 737 370 | 617 | 863 
299 | 570 | 764 339 | 570 | 795 | p=125 | 434 | 500 | 809 
298 | 574 | 767 | p=115 | 418 | 437 | 740 403 | 550 | 828 
280 | 665 | 840 399 | 460 | 744 398 | 560 | 833 
278 | 678 | 851 343 | 572 | 800 373 | 622 | 870 
273 | 714 | 882 323 | 644 | 852 372 | 625 | 872 
271 | 730 | 896 318 | 667 | 870 348 | 715 | 938 
264 | 795 | 954 | p=116 | 404 | 462 | 750 341 | 750 | 966 
260 | 840 | 995 375 | 508 | 767 | p=126 | 439 | 502 | 815 
p=106 | 336 | 477 | 707 368 | 522 | 774 425 | 522 | 821 
316 | 527 | 737 346 | 577 | 807 385 | 602 | 861 
300 | 583 | 777 345 | 580 | 809 375 | 630 | 879 
280 | 689 | 863 322 | 667 | 873 | p=127 | 441 | 508 | 822 
273 | 742 | 909 317 | 692 | 893 379 | 632 | 884 
263 | 842 | 999 299 | 812 | 995 | p=128 | 446 | 510 | 828 
p=107 | 371 | 428 | 692 | p=117 | 348 | 585 | 816 382 | 637 | 891 
319 | 532 | 744 340 | 611 | 834 381 | 640 | 893 
318 | 535 | 746 | p=118 | 352 | 587 | 821 350 | 764 | 986 
p=108 | 376 | 430 | 698 | p=119 | 413 | 476 | 770 | p=129 | 428 | 546 | 845 
321 | 540 | 753 355 | 592 | 828 400 | 602 | 873 
295 | 644 | 831 354 | 595 | 830 384 | 645 | 900 
p=109 | 378 | 436 | 705 338 | 651 | 870 | p=130 | 473 | 494 | 837 
325 | 542 | 758 | p=120 | 418 | 478 | 776 388 | 647 | 905 
297 | 654 | 842 408 | 492 | 780 375 | 690 | 935 
p=110 | 383 | 4388 | 711 391 | 520 | 791 | p=131 | 455 | 524 | 848 
328 | 547 | 765 378 | 546 | 804 391 | 652 | 912 
327 | 550 | 767 374 | 555 | 809 390 | 655 | 914 
284 | 765 | 939 364 | 580 | 824 | p=132 | 460 | 526 | 854 
p=111 | 405 | 420 | 714 357 | 600 | 837 393 | 660 | 921 
330 | 555 | 774 328 | 716 | 924 | p=133 | 473 | 518 | 858 
319 | 592 | 800 323 | 744 | 947 462 | 532 | 861 
310 | 629 | 828 | p=121 | 438 | 462 | 779 435 | 574 | 876 
294 | 717 | 900 420 | 484 | 783 418 | 608 | 893 
p=112 | 370 | 476 | 734 361 | 602 | 842 397 | 662 | 926 
366 | 483 | 737 350 | 638 | 867 390 | 684 | 941 
348 | 520 | 756 339 | 682 | 900 | p=134 | 467 | 534 | 867 
334 | 557 | 779 330 | 726 | 935 400 | 667 | 933 


Table 3.8: Nontrivial Solution Chart 


p=135 


p=136 


p=137 


p=138 


p=139 


p=140 


p=141 


p=142 


p=143 


p=144 


p=145 


p=146 


p=148 


p=149 


p=150 


p=151 


p=152 


p=153 
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Chapter 4 


Conclusion 


I was able to locate a minimum of one solution for all p values up to 153 with 
the restriction that p>2 depicted in Theorem 3.1.1. It’s also important to note that x4k, 
for non-trivial solutions, 1<k, and x<n. All triples above satisfy identity 3.3. 


We can also observe that, from the resulting values, the sum of the k and x 
values will always be larger than the corresponding n value. There will always be a so- 
lution when x>5 and k>3. Also, based on the above charts, we can observe that, for 
some p values, there exist multiple triple solutions. It is important to note that there 
are solutions in which k = 2p. Moreover, the p-value is never listed as either a k, x, or n 


value and all the solutions for k, x, and n are always larger than the p-value. 


Lastly, for the n value, there were no recognizable patterns that would aid in 
predicting future solutions aside from the noticeable fact that the values are increasing 
as the p-value increases even though there are some repetitions in values of n. How- 
ever, previous researchers discovered that we can obtain multiple non-trivial solutions by 
putting some bounds on n. For instance, an online journal published by The Mathemat- 
ical Association of America on the sum of cubes provided the n? + n + 1 4 q and n? 
+ n+ 1 ¥ 3q for some arbitrary q then we have a non-trivial solution. These equations 
were instrumental in generating the solution charts above. However, there were other 
solutions, obtained, that did not always align with the above parameters for gq and 3q; 


thus, we cannot concretely stand on a definite pattern for n. 


Appendix A 


How to Generate Faulhaber’s 


Triangle 
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How To Generate Faulhaber’s Triangle 


The C++ code provided generates Faulhaber’s Triangle values up to row 34 
when you see the phrase for (int i = 0; i < 34; i++) in the actual code. The entire code 
was created by a website, Rosetta Code. The link is provided in the bibliography page. 
An alternative code is also provided in the Java language. Any experienced coder can 
copy and paste the code into the program and the triangle will generate. Feel free to 


extend the triangle beyond row 34. 
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Language C++ 
Translation of: C(Hashtag) 
Uses C++ 17 


#include <exception> 
#include <iomanip> 
#include <iostream> 
#include <numeric> 
#include <sstream> 


#include <vector> 


class Frac { 

public: 

Frac(long n, long d) { 

if (d == 0) { 

throw new std::runtime_error("d must not be zero"); 


} 


long nn = n; 


long dd = d; 
if (mn == 0) { 
dd = 1; 


} else if (dd < 0) f 


mn = -nn; 
dd = -dd; 
} 


long g = abs(std::gcd(nn, dd)); 


if (g > 1) { 
nn /= g; 
dd /= g; 


} 


num = nn; 
denom = dd; 


as 


Frac operator-() const { 


return Frac(-num, denom) ; 


} 


Frac operatort+(const Frac& rhs) const { 


return Frac(num*rhs.denom + denom * rhs.num, rhs.denom*denom) ; 


z 


Frac operator-(const Frac& rhs) const { 


return Frac(num*rhs.denom - denom * rhs.num, rhs.denom*denom) ; 


as 


Frac operator*(const Frac& rhs) const { 


return Frac(num*rhs.num, denom*rhs.denom) ; 


7 


friend std::ostream& operator<<(std::ostreamz, const Frac); 


static Frac ZERO() { 
return Frac(0, 1); 


} 


private: 
long num; 
long denom; 


3; 
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std::ostream & operator<<(std::ostream & os, const Frac &f) { 
if (f.num == 0 || f.denom == 1) { 
return os << f.num; 


- 


std::stringstream ss; 
ss << f.num << "/" << £.denom; 


return os << ss.str(); 


} 


Frac bernoulli(int n) { 
if (n <0) { 
throw new std::runtime_error("n may not be negative or zero"); 


} 


std::vector<Frac> a; 

for (int m = 0; m <= n; m++) { 
a.push_back(Frac(1, m + 1)); 

for (int j =m; j >= 1; j--) { 

alj - 1] = (alj - 1] - alj]) * Frac(j, 1); 
5 

} 


// returns ’first’ Bernoulli number 
if (nm != 1) return a([0]; 
return -a[0]; 


} 


int binomial(int n, int k) { 
if (a <0 Jk <0 a © we) 
throw new std::runtime_error("parameters are invalid"); 


} 
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if (mn == || k == 0) return 1; 


for (int i =k +1; i <=n; i++) { 


int denom = 1; 
for (int i = 2; i <=n- k; i++) { 
denom *= i; 


hs 


return num / denom; 


us 


std::vector<Frac> faulhaberTraingle(int p) { 


std: :vector<Frac> coeffs; 


for (int i=0; i<p+41; it+) { 
coeffs.push_back (Frac: :ZERO()) ; 
: 


Frac q{.1,. ‘p+ 1) #3 

int sign = -1; 

for (int j = 0; j <= p; j++) { 

Sign *= -1; 

coeffs[p - j] = q * Frac(sign, 1) * Frac(binomial(p + 1, j), 1) * bernoulli(j); 
} 


return coeffs; 


ts 


int main() { 


using namespace std; 


for (int i = 0; i < 34; i++) f{ 

vector<Frac> coeffs = faulhaberTraingle(i) ; 

for (auto it = coeffs.begin(); it != coeffs.end(); it++) { 
cout << right << setw(5) << *it <<" "; 

} 

cout << endl; 


- 


return 0; 


As 
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How To Generate Non Trivial 


Solutions 
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How ‘To Generate Non'Trivial Solutions 


These codes were generated using the aforementioned parameters given for p>2, 
1<k, x<n, x4#k. These rules are those implemented for the resulting identity 3.3. This 
section provides the C++ codes needed to be able to generate non-trivial solutions. A 
C++ knowledgeable person is well aware of files that end in .h and those that end in 
.cpp. 

All the codes for each individual file is displayed with their titles. Be mindful 
that all .h and .cpp files must be open in order for the solutions to generate in main.cpp. 
I took the liberty of detailing the codes for each independent file on separate pages to 
avoid any confusion. Once you’ve copied and pasted the codes into C++ and made minor 
modifications when you see this phrase (move up a line) or (move up to the previous line). 
To avoid cutting off the codes at the end of each line I moved it to the next line. Make 
sure to move it back up before running the program. Once the modifications have been 
made, run the program. You should a long list of solutions as previously displayed above. 

You can extend the solutions further than what I provided by modifying the 
value listed but be mindful that there will be a considerable wait before an output is 


received. 
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Individual.cpp Code For File: 


//Individual.cpp 


#include "Individual .h" 


Individual: : Individual () 
: SumOfCubes() {} 


inline void Individual: :produceCombinations(int arr[], int data[], const int &start, 


const int &end, const int &index, const int &combinations){} (move up a line) 


//Function to set k value 

//Good! 

inline void Individual::kSetter(int arr[])f 
cout << "Please enter a value for k: "; 
cin >> this->k; 
cout << endl; 


xSetter (arr) ; 


//Function to set x value 

//Good! 

inline void Individual::xSetter(int arr[]){ 
cout << "Please enter a value for x: "; 
cin >> this->x; 
cout << endl; 


nSetter (arr) ; 


//Function to set n value 


//Good! 


inline void Individual: :nSetter(int arr[]){ 

cout << "Please enter a value for n: "; 

cin >> this->n; 

cout << endl; 

while(calc_NBoundaries() == false){ 
cout << "n does not work! Please enter a new n: "; 
cin >> this->n; 
cout << endl; 

i 

if (Boundaries() == true){ 
checkEquality() ; 

J 

else{ 


kSetter (arr) ; 


//Function to test if n is possible 
inline bool Individual: :calc_NBoundaries(){ 
int count (0); 
int q = (GetN_value()*GetN_value() + GetN_value() + 1); 
if(q == 0){ 
cout << "q cannot be equal to 0:\n"; 


return false; 


t 

if(q == 3){ 
cout << "q cannot be equal to 3:\n"; 
return false; 

7 


tf (q 4d == 0) 4 
cout << "q cannot be divisible by 3:\n"; 


return false; 


AT 
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} 
// for(unsigned int i = 2; i <q; ++i){ 


// if(q %i == 0)f 


// count++; 

ae } 

sae 

// if(count > 1){ 

// cout << "q is not prime:\n"; 
// return false; 

// } 

else{ 


return true; 


//Function to test if numbers are correct 


bool Individual: :Boundaries(){ 


long double p = (GetX_value() + GetK_value() - GetN_value()); 
if(GetK_value() != GetN_value()){ 
//Removed stipulation that x <= n for testing purposes 
if (GetX_value() != GetN_value()){ //&& GetX_value() <= GetN_value()){ 
long double xPlus_K = GetX_value() + GetK_value() ; 
if(xPlus_K > GetN_value()){ 
if (GetK_value() >= 3){ 
if (GetX_value() >= 5){ 
long double nMinusOne = GetN_value() - 1; 
if (GetK_value() != nMinusOne){ 
if ( p= 2)4 
cout << "These values work!\n"; 


return true; 


AQ 


elset{ 
cout << " p, which is the sum of (k,x,n) is not 


greater than or equal to 2.\n"; (move up a line) 


return false; 


is 
} 
else{ 
cout << "k value is equal to n value minus 1.\n"; 
return false; 
} 
} 
else{ 
cout << "x value is not greater than or equal to 5.\n"; 
return false; 
} 
t 
else{ 
cout << "k value is not greater than or equal to 3.\n"; 
return false; 
} 
A; 
else{ 
cout << "x + K is less than n value.\n"; 
return false; 
} 
} 
elsef{ 


if (GetX_value() == GetN_value()){ 
cout << "x value is equal to n value.\n"; 


return false; 
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else{ 
cout << "x value is not less than or equal to n value.\n"; 


return false; 


} 
} 
a, 
else{ 
cout << "k value is equal to n value.\n"; 
return false; 
} 


return false; 


inline bool Individual: :checkEquality(){ 


long double LHS = (Summation_KLS() + X_valueLS() + Summation_NLS()) ; 
if (LHS == Right_HSquare()){ 
cout << "The combination of this (k,x,n) - ("; 
cout << GetK_value() <<", " << GetX_value() << ", " << GetN_value() << ") 
works!\n"; (move up a line) 
return true; 
a 
else{ 
cout << "The combination of this (k,x,n) - ("; 
cout << GetK_value() <<", " << GetX_value() <<", " << GetN_value() << ") does 
not work!\n"; (move up a line) 


return false; 


Ls 


//Function to calculate the sum of k left side 
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//Good! 
inline unsigned long long int Individual: :Summation_KLS() { 
sumK_LS = 0; 
for(unsigned int j = 1; j < GetK_value(); ++j){ 
sumK_LS += j*j*j; 
5 
cout << "The summation of " << GetK_value() << " on the left-hand side is: " 
<< sumK_LS << endl; (move up a line) 


return sumK_LS; 


//Function to calculate x to thrid power 

//Good! 

inline unsigned long long int Individual: :X_valueLS(){ 
long double x_cubed(GetX_value() *GetX_value() *GetX_value()) ; 
cout << GetX_value() << " raised to the third power on the left-hand side is: " 
<< x_cubed << endl; (move up a line) 


return x_cubed; 


//Function to calculate the sum of n left-hand side 
//Good! 
inline unsigned long long int Individual: :Summation_NLS(){ 
sumN_LS = 0; 
unsigned int begin(GetK_value() + 1); 
unsigned int condition(GetN_value() + 1); 
for(unsigned int j = begin; j < condition; ++j){ 
sumN_LS += j*j*j; 
i; 
cout << "The summation of " << GetN_value() << " on the left-hand side is: " 
<< sumN_LS << endl; (move up a line) 


return sumN_LS; 
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//Function to calculate the sum of K right-hand side 
//Good! 
inline unsigned long long int Individual: :Summation_KRS() { 
sumK_RS = 0; 
for(unsigned int j = 1; j < GetK_value(); ++j){ 
sumK_RS += j; 
t 
cout << "The summation of " << GetK_value() << " on the right-hand side is: " 
<< sumK_RS << endl; (move up a line) 
cout << "The value of x on the right hand side is: " << GetX_value() << endl; 


return sumK_RS; 


//Function to calculate the sum of N right-hand Side 
//Good! 
inline unsigned long long int Individual: :Summation_NRS(){ 
sumN_RS = 0; 
unsigned int begin(GetK_value() + 1); 
unsigned int condition(GetN_value() + 1); 
for(unsigned int j = begin; j < condition; ++j){ 
sumN_RS += j; 
i, 
cout << "The summation of " << GetN_value() << " on the right-hand side is: " 
<< sumN_RS << endl; (move up a line) 


return sumN_RS; 


//Function to square right-hand side 
//Good! 


inline unsigned long long int Individual: :Right_HSquare() { 


long double RHS(Summation_KRS() + GetX_value() + Summation_NRS()); 
RHS *= RHS; 
cout << "The right-hand side summation squared is: " << RHS << endl; 


return RHS; 
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Individual.h Code For File: 


//individual.h 


/*Used to find the soltuions 
for user determined 
(k, x, n) 

*/ 


#include "SumO0fCubes.h" 


#ifndef _INDIVIDUAL_H_ 
#define _INDIVIDUAL_H_ 


class Individual : public Sum0fCubes{ 


protected: 


public: 
Individual () ; 
void kSetter(int arr[])override; 
void xSetter(int arr[]) override; 
void nSetter(int arr[]) override; 
bool calc_NBoundaries() override; 
bool Boundaries () override; 
bool checkEquality () override; 
unsigned long long int Summation_NLS() override; 
unsigned long long int Summation_KLS() override; 
unsigned long long int Summation_KRS() override; 
unsigned long long int Summation_NRS() override; 
unsigned long long int X_valueLS() override; 
unsigned long long int Right_HSquare() override; 


void produceCombinations(int arr[], int data[], const int &, 
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const int &, const int &, const int &)override; (move up a line) 


3; 


#endif 


NontTrivial.h Code For File: 


//Non-Trivial.h 


/*Used to find non-trivial 
solutions for user determined 
range for (k, x, n) 


*/ 
#include "SumOfCubes.h" 


#ifndef _NONTRIVIAL_H_ 
#define _NONTRIVIAL_H_ 


class NonTrivial : public Sum0fCubes{ 
private: 


int indexes; 


public: 

NonTrivial(); 

void produceCombinations(int arr[], int data[], const int 
const int &, const int &, const int &)override; (move up 

void kSetter(int arr[])override; 

void xSetter(int arr[]) override; 

void nSetter(int arr[]) override; 

int Amount(); 

int Amount2(); 

string Saved(int) ; 

int GetBeginRange() ; 

int GetEndRange() ; 

bool checkEquality () override; 


unsigned long long int Summation_NLS() override; 
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&, 


a line) 


3; 


#endif 


unsigned long long int Summation_KLS() override; 
unsigned long long int Summation_KRS() override; 
unsigned long long int Summation_NRS() override; 
unsigned long long int X_valueLS() override; 
unsigned long long int Right_HSquare() override; 
bool Boundaries () override; 


bool calc_NBoundaries() override; 
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Nontrivial.cpp Code For File: 


//NonTrivial.cpp 


#include "NonTrivial.h" 


NonTrivial::NonTrivial () 


:Sum0fCubes(), indexes(0) {} 


inline void NonTrivial::produceCombinations(int arr[], int data[], const int &start, 
const int &end, const int &index, const int &combinations){ (move up a line) 


//array to hold copy of range 


if (index == combinations) { 
if(data[0O] != data[1i] && data[0O] != data[2] && data[1] != data[2])f{ 
kSetter (data) ; 


// replace index with all possible 

// elements. The condition "end-i+1 >= r-index" 
// makes sure that including one element 

// at index will make a combination with 


// remaining elements at remaining positions 


for(unsigned int i = start; i <= end && (end - i+ 1 >= combinations - index); 
++i){ (move up a line) 
data[index] = arr[i]; 


produceCombinations(arr, data, i+1, end, index+1, combinations) ; 


//Function to set k value 


59 


//Good! 
inline void NonTrivial::kSetter(int arr[]) { 
this->k = arr([0]; 


xSetter (arr) ; 


inline void NonTrivial::xSetter(int arr[]){ 
this->x = arr[1]; 


nSetter (arr) ; 


inline void NonTrivial: :nSetter(int arr[]){ 
this->n = arr[2]; 
if (calc_NBoundaries() == true){ 
if (Boundaries() == true){ 
if (checkEquality() == true){ 
} 


//Function to test if numbers are correct 


inline bool NonTrivial: :Boundaries(){ 


long double p(GetX_value() + GetK_value() - GetN_value()); 
if (GetK_value() != GetN_value()){ 

//Removed stipulation that x <= n for testing purposes 

if (GetX_value() != GetN_value()){ //&& GetX_value() <= GetN_value()){ 
long double xPlus_K(GetX_value() + GetK_value()); 
if(xPlus_K > GetN_value()){ 

if (GetK_value() >= 3){ 
if (GetX_value() >= 5){ 


long double nMinusOne(GetN_value() - 1); 


if (GetK_value() != nMinusOne){ 
if( p >= 2){ 


return true; 


as 


return false; 


inline bool NonTrivial: :calc_NBoundaries(){ 
int count (0); 
int q(GetN_value()*GetN_value() + GetN_value() + 1); 
if(q == 0){ 


return false; 


i; 

if(q == 3)f 
return false; 

i; 


if(q. 403 == 0)4 


return false; 


z 
else{ 

return true; 
} 


// if(q % 7 == 0)f 


// return true; 
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// ¥4 


// elset 
// return false; 
// } 


inline bool NonTrivial: :checkEquality( { 


long double LHS(Summation_KLS() + X_valueLS() + Summation_NLS()); 
if (LHS == Right_HSquare()){ 
ListOfSolutions(GetK_value(), GetX_value(), GetN_value()); 
return true; 
ts 


return false; 


//Function to calculate the sum of k left side 
//Good! 
inline unsigned long long int NonTrivial: :Summation_KLS(){ 
sumK_LS = 0; 
for(unsigned int j = 1; j < GetK_value(); ++j){ 
sumK_LS += j*j*j; 
t 


return sumK_LS; 


//Function to calculate x to thrid power 
//Good! 
inline unsigned long long int NonTrivial::X_valueLS(){ 
long double x_cubed (GetX_value() *GetX_value() *GetX_value()) ; 


return x_cubed; 
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//Function to calculate the sum of n left-hand side 
//Good! 
inline unsigned long long int NonTrivial: :Summation_NLS(){ 
sumN_LS = 0; 
unsigned int begin (GetK_value() + 1); 
unsigned int condition (GetN_value() + 1); 
for(unsigned int j = begin; j < condition; ++j){ 
sumN_LS += j*j*j; 
} 


return sumN_LS; 


//Function to calculate the sum of K right-hand side 
//Good! 
inline unsigned long long int NonTrivial: :Summation_KRS() { 
sumK_RS = 0; 
for(unsigned int j = 1; j < GetK_value(); ++j){ 
sumK_RS += j; 
} 


return sumK_RS; 


//Function to calculate the sum of N right-hand Side 
//Good! 
inline unsigned long long int NonTrivial: :Summation_NRS() { 
sumN_RS = 0; 
unsigned int begin(GetK_value() + 1); 
unsigned int condition(GetN_value() + 1); 
for(unsigned int j = begin; j < condition; ++j){ 


sumN_RS += j; 


} 


return sumN_RS; 


//Returns the size of the vector 
int NonTrivial: :Amount(){ 
indexes = Answerlist.size(); 


return indexes; 


int NonTrivial: :Amount2(){ 
indexes = correctList.size(); 
return indexes; 
- 
//Returns the results at each index 
string NonTrivial::Saved(int index) { 


return Answerlist.at (index) ; 


int NonTrivial: :GetBeginRange(){ 


return this->beginRange ; 


int NonTrivial: :GetEndRange(){ 


return this->endRange; 


//Function to square right-hand side 
//Good! 


unsigned long long int NonTrivial: :Right_HSquare(){ 


long double RHS(Summation_KRS() + GetX_value() + Summation_NRS()); 


RHS *= RHS; 
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return RHS; 
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SumOfCubes.h Code For File: 


//SumO0fCubes.h 


//An upgraded version of this program will use pointers 


//This will improve memory management 


// #include "Buddy.h" 
#include <iostream> 
#include <array> 
#include <fstream> 
#include <string> 
#include <sstream> 


#include <vector> 


#include "SubSet.h" 


using namespace std; 


#ifndef _SUMOFCUBES_H_ 
#define _SUMOFCUBES_H_ 


class Sum0fCubes : public SubSet{ 
protected: 
unsigned short int k; 
unsigned short int x; 
unsigned short int n; 
unsigned long long int sumK_LS; 
unsigned long long int sumN_LS; 
unsigned long long int sumK_Rs; 
unsigned long long int sumN_Rs; 


int beginRange; 
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int endRange; 

int combinations; 

int SIZE_OF_ARRAY ; 

vector<string> Answerlist; 

int pCount; 

int *range; 

int *rangeDuplicate; 

short int position; 

vector <SubSet> correctList; 

vector<int> PeeCount; 
public: 

SumOfCubes () ; 

void DefineRange() ; 


virtual void produceCombinations(int arr[], int data[], const int &, 


const int &, const int &, const int &)=0; (move up a line) 
virtual void kSetter(int arr[]) = 0; 

unsigned short int GetK_value() ; 

virtual void xSetter(int arr[]) = 0; 


unsigned short int GetX_value() ; 


virtual void nSetter(int arr[]) = 0; 


unsigned short int GetN_value(); 


virtual bool Boundaries() = 0; 


virtual bool calc_NBoundaries() = 0; 


virtual bool checkEquality() = 0; 


virtual unsigned long long int Summation_NLS() 


virtual unsigned long long int Summation_KLS() 


virtual unsigned long long int Summation_KRS() 


virtual unsigned long long int Summation_NRS() 


virtual unsigned long long int X_valueLS()=0; 


virtual unsigned long long int Right_HSquare() 0; 
void ListOfSolutions(const int &, const int &, const int &); 


void Print(); 


3; 


#endif 


void GetHighestPCount (unsigned int &); 
int GetPCount (unsigned int &i); 

void EmptyVector() ; 

void deletePointerArray(); 


bool FindPattern(unsigned int &); 
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SumOfCubes.cpp Code For File: 


//Sum0fCubes.cpp 
#include "SumO0fCubes.h" 


SumO0fCubes: : Sum0fCubes () 

: SubSet(0O, 0, 0), k(O), x(O), n(0), sumK_LS(0), 
sumN_LS(0), sumK_RS(0), sumN_RS(0O), 

beginRange(0), endRange(0), combinations (3), 
SIZE_OF_ARRAY(0O), Answerlist(0), pCount(2), range(nullptr), 
rangeDuplicate(nullptr), position(0) {} 


//Function to define range of numbers 
//Good! 
void Sum0fCubes: :DefineRange() { 
cout << "Please enter the Low-end of the range: "; 
cin >> this->beginRange; 
cout << endl; 
cout << "Please enter the high-end of the range: "; 
cin >> this->endRange; 
cout << endl; 
long double size = ((this->endRange - this->beginRange) + 1); 
SIZE_OF_ARRAY = size; 
range = new int [SIZE_OF_ARRAY] ; 
for(unsigned int i = 0; i < SIZE_OF_ARRAY; ++i){ 
range[i] = this->beginRange; 
this->beginRanget+ ; 
} 
cout << endl; 


rangeDuplicate = new int[combinations] ; 


size = SIZE_OF_ARRAY - 1; 
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produceCombinations(range, rangeDuplicate, 0, size, 0, combinations) ; 


//Function to return k value 
//Good! 
unsigned short int Sum0fCubes: :GetK_value(){ 


return this->k; 


//Function to return x value 
//Good! 
unsigned short int Sum0fCubes: :GetX_value(){ 


return this->x; 


//Function to get n value 
//Good! 
unsigned short int Sum0fCubes: :GetN_value(){ 


return this->n; 


//Function to store correct answers in vector 
void Sum0fCubes: :ListOfSolutions(const int &k_value, const int &x_value, const int 


&n_value){ (move up to the previous line) 


SubSet correct(k_value, x_value, n_value); 


correctList.push_back(correct) ; 
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void Sum0fCubes: :Print(){ 


int 
int 
int 
int 


int 


iPCount ; 
jPCount ; 
location; 
HighestPCount ; 


tabSpaces(1) ; 


ostringstream stri; 


string finalString; 


for(unsigned int i = 0; i < correctList.size(); ++i){ 


GetHighestPCount (i); ; 


HighestPCount = pCount; 


pCount = 2; 


while(pCount < HighestPCount) { 


cout << "When p = " << pCount << endl; 
stri << "When p = " << pCount << endl; 
for(unsigned int i = 0; i < correctList.size(); ++i){ 
if (GetPCount(i) == pCount){ 
cout << "(" << correctList.at(i).KC << ", " << correctList.at(i) .XC 
<<", " << correctList.at(i).NC << ")" << "\t"; (move up a line) 
stri << "(" << correctList.at(i).KC << ", " << correctList.at(i) .XC 
<< "," << correctList.at(i).NC << ")" << "\t"; (move up a line) 
if(tabSpaces % 5 == 4){ 
cout << endl; 
stri << endl; 


++tabSpaces; 
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++tabSpaces; 


By 

++pCount ; 
tabSpaces = 1; 
cout << endl; 


stri << endl; 


finalString = stri.str(); 


Answerlist.push_back(finalString) ; 


inline void Sum0fCubes: :GetHighestPCount (unsigned int &i){ 


int count = correctList.at(i).XC+tcorrectList.at(i).KC - correctList.at(i). NC; 
if(count > pCount){ 


pCount = count; 


inline int Sum0fCubes::GetPCount (unsigned int &i){ 


int count = correctList.at(i).XC+correctList.at(i).KC - correctList.at(i). NC; 


return count; 


void Sum0fCubes: :EmptyVector(){ 
pCount = 0; 
correctList.clear() ; 


Answerlist.clear(); 


void Sum0fCubes: :deletePointerArray(){ 
delete [] range; 
range = nullptr; 
delete [] rangeDuplicate; 


rangeDuplicate = nullptr; 


inline bool Sum0fCubes::FindPattern(unsigned int &s){ 
int checkN = correctList.at(s).NC; 
int theorem = (checkN*checkN+checkN + 1); 
if(theorem % 7 == 0){ 


return true; 


} 
else{ 

return false; 
i, 
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Subset.h Code For File: 


//SubSet.h 


#include <vector> 


#ifndef _SUBSET_H_ 
#define _SUBSET_H_ 


class SubSet{ 
public: 
unsigned short int KC; 
unsigned short int XC; 
unsigned short int NC; 
SubSet (const int &, const int &, const int &); 
ae 
#endif} 
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Subset.cpp Code For File: 


//SubSet .cpp 
#include "SubSet.h" 


SubSet: :SubSet (const int &KCorrect, const int &XCorrect, const int &NCorrect) 


: KC(KCorrect), XC(XCorrect), NC(NCorrect) {} 
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Main.cpp Code For File: 


//main.cpp 

#include "Sum0fCubes.h" 
#include "Individual.h" 
#include "Trivial.h" 
#include "NonTrivial.h" 
#include <fstream> 
#include <cassert> 


inline void WelcomeScreen() { 


cout 


cout 


cout 


cout 


cout 


cout 


cout 


cout 


<< endl; 

<< "Welcome to Sum of Cubes! \n"; 

<< "The purpose of this program is to\n"; 

<< "find combinations (k, x, n) that make\n"; 
<< "equations on the right and left equal\n"; 
<< endl; 

<< "What would you like to do?\n"; 


<< endl; 


inline void Menu(){ 


cout 


cout 


cout 


cout 


cout 


cout 


<<" --------------------------------------------------------------------- 
<< "1. Enter (k, x, n) individually\n"; 

<< "2. Find non-trivial solutions for (k, x, n) within a given range\n"; 
<< "3. Find Trivial solutions for (k, x, n) within a given range\n"; 

<< "4. Quit\n"; 


Seo Wiese Ko eo eee ee a a ee ee eee ae oe Fe ee eee 


int main(){ 


int option(0); 

char choice; 

int fileAppendage2(0) ; 
int fileAppendage3(0) ; 
int iterator(0); 

int iterator2(0); 

int Emptyarr [3] ; 

string NonTrivialSolutions; 
string TrivialSolutions; 
string example; 
ostringstream fileName2; 


ostringstream fileName3; 


Individual Individual; 
NonTrivial NonTrivial; 


Trivial Trivial; 


ofstream outFS; 


WelcomeScreen() ; 
Menu(); 
cin >> option; 
cout << endl; 
while(option != 4){ 
if(option == 1){ 
Individual .kSetter(Emptyarr) ; 
Menu () ; 
cin >> option; 
Bs 
else 


if(option == 2){ 
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NonTrivial .DefineRange() ; 
NonTrivial.deletePointerArray() ; 
NonTrivial.Print() ; 
cout << endl; 
cout << "Would you like to save this data? (Y/N) "; 
cin >> choice; 
if(choice == ’Y’ || choice == ’y’){ 
fileName2 << "Non-Trivial Solutions" << fileAppendage2 << ".txt"; 
NonTrivialSolutions = fileName2.str(); 
outFS.open(NonTrivialSolutions) ; 
if (!outFS.is_open()){ 


cout << "could not open file!\n"; 


a: 
else{ 
iterator2 = NonTrivial.Amount2(); 
iterator = NonTrivial.Amount(); 
outFS << "The total number of solutions from " << NonTrivial. 
GetBeginRange() - NonTrivial.GetEndRange(); (move up a line) 
outFS << " to " << NonTrivial.GetEndRange() << " are: " 
<< iterator2; (move up to previous line) 
outFS << endl; 
outFS << NonTrivial.Saved(0) << endl; 
} 
t 
++fileAppendage?2; 


fileName2.str(string()); 
outFS.close() ; 
NonTrivial.EmptyVector() ; 
Menu () ; 
cin >> option; 

} 

else{ 
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Trivial .DefineRange() ; 

Trivial.deletePointerArray () ; 

Trivial.Print(); 

cout << "Would you like to save this data? (Y/N) "; 

cin >> choice; 

if(choice == ’Y’ || choice == ’y’){ 
fileName3 << "Trivial Solutions" << fileAppendage3 << ".txt"; 
TrivialSolutions = fileName3.str() ; 
outFS.open(TrivialSolutions) ; 
if (!outFS.is_open()){ 


cout << "could not open file!\n"; 


i; 
else{ 
iterator = Trivial.Amount(); 
outFS << "The total number of solutions from " 
<< Trivial.GetBeginRange(); (move up to previous line) 
outFS << " to " << Trivial.GetEndRange() << "are: " << iterator; 
outFS << endl; 
for(unsigned int i = 0; i < iterator; ++i){ 
outFS << Trivial.Saved(i) << endl; 
} 
} 
} 
++fileAppendage3; 


fileName3.str(string()) ; 
outFS.close() ; 
Trivial.EmptyVector() ; 
Menu() ; 


cin >> option; 


return 0; 
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Appendix C 


How To Generate Trivial 


Solutions 
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How ‘To Generate Trivial Solutions 


This wasn’t a focus in the research because it’s relatively easier to generate than non- 
trivial solutions. However, for the overzealous mathematicians who would like to see what 
trivial solutions look like. I provided the C++ code below. Feel free to make comparisons 


between Non-trivial and Trivial solutions. 
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Trivial.cpp Code For File: 


//Trivial.cpp 


#include "Trivial.h" 


//constructor 
Trivial: :TrivialQ 


: SumOfCubes(), indexes(0) {} 


inline void Trivial: :produceCombinations(int arr[], int data[], const int &start, 


const int &end, const int &index, const int &combinations){ (move up a line) 


if (index == combinations) { 
if(data[0] == data[1i] || data[0] == data[2] || data[1] == data[2]){ 
cout << endl; 


kSetter (data) ; 


// replace index with all possible 

// elements. The condition "end-i+1 >= r-index" 
// makes sure that including one element 

// at index will make a combination with 


// remaining elements at remaining positions 


for(unsigned int i = start; i <= end && (end - i+ 1 >= combinations - index); 
++i){ (move up a line) 
data[index] = arr[i]; 


produceCombinations(arr, data, i+1, end, index+1, combinations) ; 


//Function to set k value 

//Good! 

inline void Trivial::kSetter(int arr[]){ 
this->k = arr[0]; 


xSetter (arr) ; 


inline void Trivial::xSetter(int arr[]){ 
this->x = arr[1]; 


nSetter (arr) ; 


inline void Trivial: :nSetter(int arr[]){ 
this->n = arr[2]; 
if (calc_NBoundaries() == true){ 
if (Boundaries() == true){ 
if (checkEquality() == true){ 
iy 


//Function to test if numbers are correct 


inline bool Trivial: :Boundaries(){ 


long double p = (GetX_value() + GetK_value() - GetN_value()); 
long double xPlus_K = GetX_value() + GetK_value() ; 
if(xPlus_K > GetN_value()){ 
if (GetK_value() >= 3){ 
if (GetX_value() >= 5){ 
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long double nMinusOne = GetN_value() - 1; 


if (GetK_value() != nMinusOne){ 
if( p >= 2){ 


return true; 


} 


return false; 


inline bool Trivial: :calc_NBoundaries(){ 
int count = 0; 
int q = (GetN_value()*GetN_value() + GetN_value() + 1); 
if(q == 0)f 


return false; 


J 

if(q == 3){ 
return false; 

J 


if(q % 3 == o)f 


return false; 


F 
else{ 

return true; 
5 


inline bool Trivial: :checkEquality(){ 
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long double LHS = (Summation_KLS() + X_valueLS() + Summation_NLS()) ; 


if (LHS == Right_HSquare()){ 
ListOfSolutions(GetK_value(), GetX_value(), GetN_value()); 
return true; 

J 


return false; 


//Function to calculate the sum of k left side 
//Good! 
inline unsigned long long int Trivial: :Summation_KLS(){ 
sumK_LS = 0; 
for(unsigned int j = 1; j < GetK_value(); ++j){ 
sumK_LS += j*j*j; 
J 


return sumK_LS; 


//Function to calculate x to thrid power 
//Good! 
inline unsigned long long int Trivial: :X_valueLS(){ 
long double x_cubed = (GetX_value()*GetX_value()*GetX_value()) ; 


return x_cubed; 


//Function to calculate the sum of n left-hand side 
//Good! 
inline unsigned long long int Trivial: :Summation_NLSQ{ 
sumN_LS = 0; 
unsigned int begin = GetK_value() + 1; 
unsigned int condition = GetN_value() + 1; 


for(unsigned int j = begin; j < condition; ++j){ 
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sumN_LS += j*j*j; 
J 


return sumN_LS; 


//Function to calculate the sum of K right-hand side 
//Good! 
inline unsigned long long int Trivial: :Summation_KRS(){ 
sumK_RS = 0; 
for(unsigned int j = 1; j < GetK_value(); ++j){ 
sumK_RS += j; 
i, 


return sumK_RS; 


//Function to calculate the sum of N right-hand Side 
//Good! 
inline unsigned long long int Trivial: :Summation_NRS(){ 
sumN_RS = 0; 
unsigned int begin = GetK_value() + 1; 
unsigned int condition = GetN_value() + 1; 
for(unsigned int j = begin; j < condition; ++j){ 
sumN_RS += j; 
i, 


return sumN_RS; 


//Function to square right-hand side 
//Good! 


inline unsigned long long int Trivial: :Right_HSquare(){ 
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long double RHS = (Summation_KRS() + GetX_value() + Summation_NRS()); 


RHS *= RHS; 


return RHS; 


//Returns the size of the vector 
int Trivial: :Amount(){ 
indexes = Answerlist.size(); 


return indexes; 


string Trivial: :Saved(int index){ 


return Answerlist.at (index) ; 


int Trivial: :GetBeginRange() { 
return this->beginRange ; 

Bg 

int Trivial: :GetEndRange(){ 


return this->endRange; 
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Trivial.h Code For File: 


//Sum0fCubes2.h 


#include "SumO0fCubes.h" 


#ifndef _TRIVIAL_ 
#define _TRIVIAL_ 


class Trivial 


private: 


H_ 
H_ 


int indexes; 


public: 


Trivial(); 


: public Sum0fCubes{ 


void produceCombinations(int arr[], int datal[], 


const int &, 
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const int &, const int &, const int &)override; (move up a line) 


int Amount(); 


string Saved(int) ; 


void kSetter(int arr[]) override; 


void xSetter(int arr[]) override; 


void nSetter(int arr[])override; 


int GetBeginRange() ; 


int GetEndRange() ; 


bool checkEquality () override; 


unsigned 
unsigned 
unsigned 
unsigned 
unsigned 


unsigned 


long long 
long long 
long long 
long long 
long long 
long long 


int Summation_NLS() override; 
int Summation_KLS() override; 
int Summation_KRS() override; 
int Summation_NRS() override; 
int X_valueLS() override; 


int Right_HSquare() override; 


bool Boundaries () override; 


bool calc_NBoundaries() override; 


a 
#endif 
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